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An improved method is proposed for the extraction of the symmetry energy coefficient relative to the tem- 
perature, dsym/T, for the heavy-ion reactions near the Fermi energy region, based on the modified Fisher Model 
(MFM). This method is applied to the primary fragments of the Anti-symmetrized Molecular Dynamics (AMD) 
simulations for the reactions of *°Ca + *°Ca at 35 MeV/nucleon. The density and the temperature at the fragment 
formation stage are extracted using a self-consistent method. 


Keywords: Intermediate energy heavy ion reactions, Symmetry energy, Density, Temperature, Modified Fisher model, Self- 


consistent method 


DOI: 10.13538/j.1001-8042/nst.26.S20508 


I. INTRODUCTION 


The symmetry energy term in the nuclear equation of state 
(EOS) intimately relates to a wealth of the dynamical pro- 
cess of nuclear reactions, the structure character of nuclei and 
astrophysical phenomena [1]. The investigations of the sym- 
metry energy, especially focusing on its density dependence, 
have been conducted through many observables such as iso- 
topic ratio [2], isospin diffusion [3], neutron-proton emission 
ratio [4, 5], giant monopole resonance [6], pygmy dipole res- 
onance [7], giant dipole resonance [8], collective flow [9] and 
isoscaling [10-12]. In our recent works, we utilized the iso- 
topic yields to extract the symmetry energy coefficient rela- 
tive to the temperature, dsym/T, in the framework of the nu- 
clear phase transition theory [13-17]. Along this scenario, 
Asym/T was extracted, using m-scaling [18] and isobaric yield 
ratios [19]. Ono et al. also independently introduced a gen- 
eralized free energy, K(N, Z), and extracted the dsym/T values 
from their quadratic distributions [20]. In this method, for a 
given Z, all contributions from the volume, surface, Coulomb 
and pairing terms in the free energy are squeezed into an N 
proportional term and a constant, &(Z)N + 7(Z), as a course 
approximation. 

While temperature is one of the key variables in charac- 
terizing nuclear reactions, it is very difficult to determine 
the temperature of hot nuclear matter in a dynamical pro- 
cess. Several nuclear thermometers have been proposed [21]. 
These include the slope of energy spectra [22, 23], momen- 
tum fluctuations [24, 25], double isotope yield ratios [26] and 
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excited state distributions [27]. However they may not be 
generally applicable in all circumstances and even for a giv- 
en system the extracted temperature values from these ther- 
mometers may be quite different from each other [21]. In 
our recent work, temperature is evaluated from a (un — Hp)/T 
analysis, but no density determination was possible [28]. In 
our another recent work, the isotopic yield ratio method was 
applied to extract dsym/T values from the experimentally re- 
constructed primary fragment yields. These ratios were com- 
pared to those calculated from AMD primary generated frag- 
ment yields obtained using Gogny interactions with different 
density dependencies of the symmetry energy [29, 30]. In the 
analysis, we found that the extracted dsym/T values change 
according to the interactions used. From the dependence on 
interaction, the density, symmetry energy and temperature at 
the time of fragment formation were determined in a self- 
consistent manner. 


In this article, an improved method is used to extract the 
dsym/T values. In the improved method, all available iso- 
tope yields are employed. The improvement is made pos- 
sible by taking into account the mass dependent tempera- 
ture [31] in the free energy in an iterative manner. Using 
the self-consistent manner, the density and the temperature 
at the time of fragment formation are carefully determined 
from the obtained asym/T values from different interactions 
with different density dependencies of the symmetry ener- 
gy term, i.e., the standard Gogny interaction which has an 
asymptotic soft symmetry energy (g0), the Gogny interac- 
tion with an asymptotic stiff symmetry energy (g0AS) and the 
Gogny interaction with an asymptotic super-stiff symmetry 
energy (gOASS) [12, 32]. Our present analyses are conducted 
in the framework of AMD. There are three major reasons to 
use AMD as the event generator for this work. One is its ca- 
pability to reproduce the experimental isotope yields. AMD 
results, such as multiplicity, angular distribution and energy 
spectra, have often been compared with those from the exper- 
imental data for intermediate energy heavy ion collisions and 
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reproduce them reasonably well [32-39]. In one of our recent 
works in Ref.[30], the yields of the experimentally recon- 
structed primary hot isotopes are well reproduced by those of 
the AMD simulations. Second is to eliminate the secondary 
cooling effect. As shown in Ref.[19, 30], the sequential decay 
of the primary hot isotopes significantly alters the yield dis- 
tribution and distorts the information inherent in the primary 
hot fragment yields. Third is to simplify the initial conditions 
using zero impact parameter to eliminate the effects of trans- 
verse flow and neck emission among others. Thus the AMD 
events are generated for central collisions (b = 0 fm) of Ca 
+ “Ca at 35 MeV/nucleon. 


I. IMPROVED MFM MODEL AND EXTRACTION OF 
Asym/ T 


In the framework of MFM, the yield of an isotope with 
mass A and J = N — Z (N neutrons and Z protons) produced 
in a multifragmentation reaction, can be given as [14, 16- 
19, 40] 


WU, A) + UnaN + pZ 


T 
(1) 
+N In) + ZÉ). 


YU, A) =Yo- A‘ exp[ 


Using the generalized Weiszdcker-Bethe semiclassical 
mass formula [41, 42], WU, A) can be approximated as 


ZZ — 1) 
WU, A) =ayA = asA23 = Go Als 
(N -Z) 6 
— asym A ~ Sp Ai? (2) 
= (ee 
ô= 7 ; 


In Eq.(1), A™ and Nin(N/A) + Zin(Z/A) originate from the 
increases of the entropy and the mixing entropy at the time of 
the fragment formation, respectively. Hn (Hp) is the neutron 
(proton) chemical potential. t is the critical exponent. In this 
work, the value of t = 2.3 is adopted from the previous stud- 
ies [40]. In general coefficients, ay, 4s, sym, 4p and the chem- 
ical potentials are temperature and density dependent. In this 
formulation a constant volume process at equilibrium is as- 
sumed in the free energy and therefore the term “symmetry 
energy” is used throughout this work along Ref. [43]. If one 
assumes a constant pressure at the equilibrium process [44], 
the term “symmetry enthalpy” should be used. Experimen- 
tally, whether the equilibrium process takes place at constant 
pressure or volume can not be determined, and thus we use 
“symmetry energy” throughout the paper, keeping in mind 
the ambiguity [43]. 

In the previous analyses [19, 30, 45, 46], the temperature 
in Eq.(1) was assumed to be identical to the temperature of 
the fragmenting source and treated as a constant for all iso- 
topes. However as seen below, this temperature turns out to 
be (fragment) mass dependent. This mass dependence on the 
temperature was not recognized in the previous analyses, just 
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because the mass dependence was masked by the larger er- 
ror bars. In this improved method, the error bars become s- 
mall and the mass dependence becomes evident. In order to 
take into account the mass dependence of the temperature in 
Eq.(1), the temperature T is replaced by an apparent temper- 
ature T(A). We attribute this mass dependence to the system 
size effect as discussed in Section IV. In the improved MFM 
formulation, therefore, this system size effect is empirically 
realized by reducing the apparent temperature as A increases 
as T(A) = To(1—-kA). To is the temperature of the fragmenting 
source and k is a constant quantifying the mass dependence. 

In order to study the density, temperature and symmetry en- 
ergy in the fragmenting source, the improved MFM of Eq.(1) 
is utilized to extract the dsym/T value from the available iso- 
tope yields. Since the dsym/To value in Eqs.(1) and (2) de- 
pends on 5 parameters, ay, ds, Ac, dp and Ap (defined by 
AH = Mn — Hp), the optimization process of these parame- 
ters is divided into the following three steps to minimize the 
ambiguity for each parameter. For a given k value 


1. Optimize Au/To and ac/To values from mirror isobars 
and fix these parameter values. 


2. Optimize a,/To, as/To and a,/To values from N = Z 
isotopes. 


3. Using extracted parameters in step (1) and step (2), 
Asym/To values are extracted from all available iso- 
topes. Comparing the extracted asym/To values with 
three different interactions, the density of the fragment- 
ing source is extracted. Using this density, the value of 
the symmetry energy coefficient, asym, for each interac- 
tion is determined. The temperature is then calculated 
as the ratio Of asym to dsym/To. 


It is expectable that if the k value is properly selected which 
means the mass dependence is well considered, a constant 
To is obtained. Since the k value is small as seen below, we 
perform the optimization of the parameter k in an iterative 
manner in the following analysis, that is, in the first round 
k = kı = Ois set in T(A) = To(1 — kA) and calculate the 
temperature as a function of A, using steps (1)-(3). From this 
plot a new k value, k = k’, is extracted from the slope. In the 
second round, k = k) = kı + $k; is used for the steps (1)-(3) 
and anew k value, k = ks, is extracted. If ky is O within a given 
error range, the iteration stops and the kz value is fixed as the 
mass dependent parameter of the apparent temperature and 
the To value is determined. Otherwise the iteration continues. 

The details of steps (1)-(3) are first described below for a 
given k value. In the step (1), following Ref. [19], the isotope 
yield ratio between isobars with J + 2 and J, RU + 2, I, A), is 
utilized, which is 


R (1+2,1,A) = YU +2,A)/YU,A) 
= exp{[Hn — Mp + 2a-(Z — 1)/A'? 

— 4asyn(I + 1)/A - 6(N + 1,Z- 1) 

— 5(N, Z)]/[To(1 — kA)] + AU + 2,1, A)}, (3) 


where Y(/, A) is the yield of isotopes with J and A, and A(Z + 
2,1,A) = Smix(+2,A)—S mixU, A). When the above equation 


$20508-2 


SYMMETRY ENERGY EXTRACTION FROM... 


i 

= 0.5) - : 4 
oot P z ] 
X e | 
= 
ate o fat 7 

L l=3 | 
NBL | E EE] 
0-3 10. 2 3 40 


Fig. 1. (Color online) In[R(I, —I, A)]/Z versus A for I = 1 and J = 3 from 
the events with the g0 interaction. The curve is the fit result of Eq.(4) for 
k = 0. The extracted values of Au/To and a,/To are given in the third and 
fifth columns of Table 1. 


is applied for a pair of mirror nuclei of odd mass isotopes with 
I = -I and J, the symmetry energy term, pairing term and 
mixing entropy terms drop out and the following equation is 
obtained. 


In[RU, -1,A)]/T = [Au + ac(A — 1)/A'*]/[To(1 — kA]. (4) 


For all available mirror isobars, Au/To and a./To are opti- 
mized in Eq.(4). The In[RU, —/, A)]/J values and the fit result 
for k = 0 is shown in Fig.1. In the step (2) we apply Eq.(1) to 
the isotopes with N = Z with the extracted Au/To and ac/To 
values in the step (1). For the N = Z = A/2 isotopes, the ratio 
of the free energy relative to the temperature can be calculat- 
ed from Eq.(1) and Eq.(2) without the symmetry energy term 
as 


_F(A/2,A/2) __ F(A/2,A/2) 


“(1 -= kA) 
To T(A) 
Y(A/2, A/2)A™ 
s OAA 
1 (5) 
-Dy _ As 42/3 _ 4c A(A= 2) 
To To To 4A1/3 
ap ô 1 
-2L A= eine 
n a tACl ~kA)IN(5), 
where ad, = a, + $ (Mn + up). The value of Inf ELADA] on 


the right of the second equation can be calculated from the 
simulated or experimental values when the t value is fixed. 
Non-zero values show the deviation of the mass distribution 
of N = Z isotopes from the power law distribution determined 
by the critical exponent [40]. When other t values are used, 
the parameter values change accordingly. In order to elimi- 
nate the constant Yo, all isotope yields are normalized to the 
yield of RC [18, 19, 40]. For the first round with k = 0, 
the renormalized values of Aare) from the AMD events 
with the g0 interaction are plotted as a function of the iso- 
tope mass A using solid points in Fig. 2(a). The values of 
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a,/To, as/To and a,/To are used as free parameters to fit the 
given fag) values, employing Eq.(5). A typical search 
result is shown by open circles in Fig. 2(a) for the case of 
the g0 interaction. Similar quality results are obtained for 
the events generated using the g0AS and gOASS interactions. 
One should note that the value of a,/T makes a small con- 
tribution and that the contribution is evident as staggering in 
the —F(A/2,A/2)/To vs. A plot. Therefore the essential free 
parameters in this step are a, /Tp and a;/To. The extracted pa- 
rameter values are summarized in Table | for the first round 
(k = 0) and the final round (k = 0.007). 


TABLE 1. a/To and Au/To for the first round (k = 0.0) and the final 
round (k = 0.007). 


k interaction 4,/Ty a;/Ty  a¢/To ap/To = Ap/To* 

0.0 g0 1.77 2.74 0.1040 0.427 -0.254 
gOAS 1.76 2.66 0.1080 0.661 -0.272 
gOASS 1.77 2.72 0.1030 0.778 -0.252 

0.007 g0 1.55 2.29 0.0805 0.416 -0.189 
gOAS 1.53 2.20 0.0830 0.651 -0.202 
gOASS 1.55 2.29 0.0800 0.766 -0.186 


a Au /To values are taken from the step (1). 


In step (3), Eq.(1) is applied to the yields of all isotopes 
with N = Z and N # Z. From Eq.(1), dsym/To and Au/To 
values can be related to the modified free energy, aN 


To 


AF(N,Z) _dsym(N-Z)* Au 


N-Z), 6 
To To A IT, ) o 


where 
F(N.Z) 
To 


AEW) is the free energy relative to the temperature, 


, subtracted by the calculated contributions of the vol- 
ume, surface, Coulomb and pairing terms, using the parame- 
ters in Table 1. Resultant arwa values are shown by symbols 
in Fig. 2(b). They exhibit quadratic relationships with mini- 
mum values close to zero. The minimum values are at or near 
N = Z isotopes and therefore reflect approximately the differ- 
ence between the data and fitting points in Fig. 2(a). In this 
step, the dsym/To and the Au/To values are optimized. Since 
the Au/To values are extracted from the step(1), the optimiza- 
tion is made around the values in the fifth column of Table 1 
in a small margin. The dsym/To values are extracted from the 
quadratic curvature of the isotope distribution for each giv- 
en Z and plotted in Fig. 2(c) separately for the g0, gOAS and 
gOASS interactions. As one can see for the first round with 
k = 0, the extracted dsym/To values increase as Z increases in 
all cases, and they more or less parallel each other. 


II. SELF-CONSISTENT DETERMINATION OF DENSITY 
AND TEMPERATURE 


In order to determine the density and temperature at the 
time of the fragment formation, the parallel behavior of the 
observed dsym/T values in Fig. 2(c) is utilized. As sug- 
gested in Ref. [12], the observed differences are attributed to 


$20508-3 


LIU Xing-Quan et al. 


E 
r t Te: J 
a T- i J 

E e° J 
wf te ' go ] 
= ot ABs ge AARTA J 
= H 24° ° 
= B a- Bi 
w f o ] 

p EET eee ESETERE 
= i 
N | 
= 7 
— 
ty J 
< J 

10; — 1 

a O 3 

E (EE + 
© [ TETI aa ] 
E 6r ms ee gotten” s J 

eun 
i =e sk ro its J 
4- 7 aay r J 
D 5 10 15 20 


Fig. 2. (Color online) (a) Calculated ratio of free energy to To for N = Z 
isotopes from the AMD events with gO (solid points). Circles represent the 
fit using Eq.(5). (b) Calculated ARN 2) values for the g0 interaction and 
quadratic fits using Eq.(6) for Z = 2 to 18. The same symbols are used for 
isotopes with a given Z. (c) Extracted dsym/To values from (b) for g0 (dots), 
gOAS (squares) and gOASS (triangles). Errors are from the quadratic fits. All 
results are from the first round (k = 0). 


the difference of the symmetry energy at the density at the 
time of the fragment formation. The ratios between g0/g0AS 
and g0/g0ASS of asym/To for the first round are shown in 
Fig. 3(a). The ratios show flat distributions as a function of Z 
for both cases. The extracted average ratio values are shown 
by lines in the figure and the values are given in the first col- 
umn of Table 2. In Fig. 3(b) the symmetry energy coefficient 
is plotted as a function of the density for the three interactions 
used in the calculations and in Fig. 3(c) their ratios, Rsym = 
Asym(g0)/dsym(gOAS ) and Rsym = Asym(g0)/dsym(gOASS ), are 
plotted. Using the ratio values determined from Fig. 3(a) and 
the density dependence of the Rsym values in Fig. 3(c), the im- 
plied densities of the fragmenting sources are indicated by the 
shaded vertical areas shown in Fig. 3(c). The extracted den- 
sity values for each case are given in the second column of 
Table 2. Assuming that the nucleon density should be same 
for the three different interactions used, the nucleon densi- 
ty of the fragmenting source is determined from the overlap 
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Fig. 3. (Color online) (a) The ratios of the asym/T values shown in Fig. 2(c), 
dots for g0/g0AS and circles for g0/g0ASS. (b) Symmetry energy coefficient 
vs density used in the AMD simulations. In the plot, solid, dashed and dotted 
lines represent the g0, gOAS and gOASS interactions, respectively. (c) The 
ratio of the symmetry energy coefficient in (b). The shaded horizontal lines 
indicate the ratios extracted in (a) and the vertical shaded area shows the 
density region corresponding to these ratios. Two different shadings are used 
for the two ratio values. All results are from the first round. 


of the extracted values. This assumption is reasonable for the 
central collisions because the nucleon density is mainly deter- 
mined by the stiffness of the EOS and not by the density de- 
pendence of the symmetry energy term. From the overlapped 
density area in Figs. 3(c), p/po = 0.67 + 0.02 is extracted as 
the density at the time of the fragment formation. Using this 
density value, the corresponding symmetry energy values at 
that density are extracted for the three different interactions 
from Fig. 3(b). They are given in the third column of Table 2. 


Once the symmetry energy value is determined for a giv- 
en interaction, the temperature, To, can be calculated as To = 
asym / (Asym/To). The extracted To values are shown as a func- 
tion of Z by open symbols for the first round in Fig. 4. The 
larger errors of Ty, comparing to those in Fig. 2(c), originate 
from the errors of asym and dsym/To extracted for each inter- 
action which are shown in the third column of Table 2 and 
Fig. 2(c), respectively. The temperature values extracted from 
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TABLE 2. Symmetry energy and p/p from k = 0.0 to k = 0.007 


k interaction Rsym p/po asym (MeV) 
0.0 go 26.6+0.3 
g0/g0AS 1.14+0.02 0.71+0.05 
g0AS 23.741.3 
g0/g0ASS 1.41+0.03 0.66+0.02 
gOASS 18.7+0.7 
0.007 go 26.8+0.2 
g0/g0AS 1.14+0.02 0.72+0.05 
g0AS 24.041.3 
g0/g0ASS 1.40+0.03 0.66+0.02 
gOASS 18.7+0.7 
8— m r 
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T, (MeV) 


Fig. 4. (Color online) Extracted To values as a function of Z. Open symbols 
are for the first round (k = 0) and closed symbols are for the final round (k = 
0.007). Different symbols represent results with g0(circles), g0AS (squares) 
and g0AS S (triangles) interactions. 


the three different interactions agree with each other very well 
and show a monotonic decrease as Z increases from ~ 5 MeV 
at Z = 4to ~ 3 MeV at Z = 18. From this slope the extract- 
ed temperature as a function of A, Ty = 5.5(1 — 0.012A), is 
determined for the first round, assuming A ~ 2Z. 

The iteration is repeated four times in this work. The same 
plots as Fig.2, but with the k value for the fourth (final) round, 
k = 0.007, are shown in Fig.5 and the extracted parameters 
are also given in Table 1. A very similar quality of results to 
those of the first round with k = 0 were obtained, even though 
the optimized parameter values are quite different between 
those of the first round (k = 0) and of the fourth round (k = 
0.007). The extracted dsym/To values parallel each other and 
show a rather flat distribution as a function of Z for Z up to 
15 in Fig. 5(c). As seen in Fig. 4, in which the extracted 
To values are shown by closed symbols as a function of Z 
for the fourth round with k = 0.007, the extracted To values 
are consistent with 5.5 MeV within the error bars. Since To 
values show a flat distribution as a function of Z, the iteration 
is stopped at this round and Tp = (5.5 + 0.2) MeV is taken as 
the temperature of the emitting source. 
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Fig. 5. (Color online) Same plots as Fig.2, but for the final round (k = 


0.007). See also the figure caption of Fig.2. 


The extracted parameter values of Rsym, density and sym- 
metry energy for the fourth round are very similar to those of 
the first round as shown in Table 2. The values and errors of 
these parameters are essentially determined by the ratios and 
their errors of the dsym/To values between different interac- 
tions as discussed in Fig.3(a). These ratio values are stable be- 
tween the first and fourth rounds, even though the optimized 
parameter values in Table 1 are quite different between these 
rounds. The extracted value of the mass dependent factor k 
has some errors, but the above fact ensures that the extracted 
density, temperature and symmetry values and their errors in 
Table 2 are rather stable independent of the choice of the k 
values within its error bar, when the parameters in Table 1 are 
optimized for the given k value. 


The decreasing trend of the temperature as A increases is 
often observed in heavy ion collisions and normally attributed 
to variations in impact parameter [47]. The heavier fragments 
tend to be produced in more peripheral collisions and there- 
fore show a lower temperature. However in this study all 
events analyzed are generated in the same class of events, 
central collisions with b = 0 fm. Therefore we attribute the 
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decreasing trend may originate to different fragment forma- 
tion processes rather than to the centrality of the events. 


IV. SUMMARY 


An improved method is proposed for the extraction of the 
ratio of the symmetry energy coefficient relative to the tem- 
perature, dsym/T(A), taking into account the mass dependence 
of the apparent temperature, based on the MFM model. This 
method is applied for the central collisions of the AMD events 
generated for *°Ca + *°Ca at 35 MeV/nucleon. The Gogny in- 
teractions, g0, gOAS and gOASS, with three different density 
dependencies of the symmetry energy are employed. As a 
function of IMF charge Z, the ratios of the extracted dgym/To 
values from different interactions are essentially constant and 
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